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ABSTRACT 

We present the analysis of six XMM-Newton observations of the Anomalous X-ray Pulsar CXOU J010043. 1- 
72 1 1 34, the magnetar candidate characterized by the lowest interstellar absorption. In contrast with all the other 
magnetar candidates, its X-ray spectrum cannot be fit by an absorbed power-law plus blackbody model. The 
sum of two (absorbed) blackbody components with kT\ = 0.30 ±0.02 keV and kT2 = 0.7 ±0.1 keV gives 
an acceptable fit, and the radii of the corresponding blackbody emission regions are R^gi 



R' 



BB2 



12.1+2-1 km and 

1 .7+Q5 km. The former value is consistent with emission from a large fraction of a neutron star surface 



and, given the well known distance of CXOU JO 10043. 1-721 134, that is located in the Small Magellanic Cloud, 
it provides the most constraining lower limit to a magnetar radius ever obtained. A more physical model, where 
resonant cyclotron scattering in the magnetar magnetosphere is taken into account, has also been successfully 
applied to this source. 

Subject headings: pulsars: individual (CXOU JO 10043. 1-72 11 34)- stars: neutron 



1. INTRODUCTION 

The Anomalous X-ray Pulsars (AXPs; see Mereghetti 
2008, for a recent review) were initially identified as a sub- 
class of accreting X-ray pulsars. Their much softer X-ray 
spectrum and the lack of a bright optical counterpart distin- 
guished them from the vast majority of X-ray pulsars, which 
are neutron stars accreting in high mass X-ray binary systems. 
AXPs have rotation periods of several seconds and show a 
secular spin-down on timescales of 10^-10^ years, but their 
rotational energy loss is smaller than their X-ray luminosity, 
excluding the possibility that they are rotation-powered, like 
radio pulsars. It is generally believed that the AXPs, as well as 
another small class of high-energy sources with similar prop- 
erties, the Soft Gamma-ray Repeaters (SGRs), are magnetars, 
i.e. neutron stars powered by their extremely high magnetic 
field ('^lO'^ G; Duncan & Thompson 1992; Thompson & 
Duncan 1996). 

The soft X-ray (1-10 keV) spectra of magnetars cannot 
be adequately fit with single component models whenever 
data with good count statistics are available. Successful fits 
are instead obtained by a two component model consisting 
of a steep power-law (photon index ~3-4) and a blackbody 
(kT ~ 0.5 keV). Some attempts have been done, also based on 
phase-resolved spectroscopy, to attribute the two components 
to physically distinct processes (e.g. Tiengo et al. 2005), but 
no particularly compelling interpretations could be obtained. 
One problem of this model is that it tends to give best-fit val- 
ues of the interstellar absorption higher than those indepen- 
dently estimated in other ways (e.g. Durant & van Kerkwijk 
2006). Another problem is that the power-law component 
cannot be extrapolated at lower energies without exceeding 
the flux of the near infrared (NIR) and optical counterparts 
(e.g. Hulleman et al. 2004). Drastic, and possibly unphysical. 
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cut-offs in the power-law component are required to match 
the low optical/NIR fluxes. 

In some AXPs good spectral fits are obtained with the sum 
of two blackbody components with different temperatures. 
Since this model does not suffer of the problems described 
above, it is usually preferred to the power-law plus blackbody 
model (e.g. Halpern & Gotthelf 2005). However, also this 
model is only phenomenological and it is inadequate to repre- 
sent the non-thermal phenomena that are expected to occur in 
the highly magnetized magnetosphere of magnetars (e.g. Lyu- 
tikov & Gavriil 2006). More physical models of the X-ray 
spectra, including the effects of the strong magnetic field and 
charged currents, have recently been developed and success- 
fully applied to a sample of magnetar candidates (Fernandez 
& Thompson 2007; Guver et al. 2007; Rea et al. 2008). 

From a purely observational point of view, it has not been 
possible to discriminate between the different models repro- 
ducing the magnetar X-ray spectra. This is mainly due to the 
low sensitivity of hard X-ray detectors above ^10 keV and to 
the large uncertainties in the fits introduced by the high inter- 
stellar absorption, that severely suppresses the flux below 1 
keV. Being young neutron stars born from massive progeni- 
tors, all the Galactic magnetars are located in highly absorbed 
regions of the Galactic plane. All of them have column den- 
sities Nh ranging from ^5 x 10^' to ^10^^ cm"^. The two 
known magnetars in the Magellanic Clouds, being consider- 
ably less absorbed, offer the possibility to better constrain the 
spectra in the low energy range. The study of SGR 0526-66, 
located in the Large Magellanic Cloud, is complicated by the 
presence of the surrounding supernova remnant N49, which 
is particularly bright in soft X-rays (Kulkarni et al. 2003). 
Here we concentrate therefore on the spectral properties of 
CXOU JO 10043. 1-72 11 34 (Lamb et al. 2002, 2003; McGarry 
et al. 2005), the only known AXP in the Small Magellanic 
Cloud (SMC). 
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TABLE 1 

Log of the XMM-Newton Observations of 
CXOU JO 1 0043 . 1-72 1 1 34 . The Pulse Periods and 
Corresponding Icr Errors Are Also Reported. 



Obs. 


Date 


PN exposure 


MOS exposure 


Period 






(ks) 


(ks) 


(s) 


A 


17 Oct 2000 


14 


20 


8.019(1) 


B 


20 Nov 2001 


22 


27 


8.0193(9) 


C 


27 Mar 2005 




24 


8.0215(9) 


D 


27 Nov 2005 


14 


17 


8.022(1) 


E 


29 Nov 2005 


13 


16 


8.022(1) 


F 


11 Die 2005 


9 


16 


8.022(2) 



^The PN data were not considered, since the PN instrument was operated 
with the filter wheel in closed position. 



2. OBSERVATIONS AND DATA ANALYSIS 

The field containing CXOU JO 10043. 1-72 1134 was ob- 
served six times'' with the XMM-Newton satellite (see Table 
1). Here we report the analysis of the data collected by the 
EPIC instrument, which is composed by one PN (Striider et al. 
2001) and two MOS X-ray cameras (Turner et al. 2001). 

CXOU JO 10043. 1-72 11 34 was not the main target of the 
observations, but, being at an off-axis angle of ~6', it was 
always well inside the field of view of the EPIC instrument 
(~15' radius). All the observations were performed with 
the medium optical blocking filter and in Full Frame mode 
(time resolution of 73 ms and 2.6 s for the PN and MOS, re- 
spectively), except for the first PN observation, done in Ex- 
tended FuU Frame mode (time resolution of 200 ms). All the 
data were processed using the XMM-Newton Science Analy- 
sis Software (SAS version 7.1.0) and the calibration files re- 
leased in August 2007. The standard pattern selection criteria 
(patterns 0^ for PN and 0-12 for MOS) were adopted. 

Source spectra were extracted for each observation from 
circular regions with 25" radius. The background spec- 
tra were extracted from a 94"x72" box centered at RA = 
Oli'OO^Se^S, Dec. = -72°11'33" and rotated such that it 
intercepts no CCD gaps in any PN image. Response matrices 
and ancillary files for each spectrum were produced using the 
SAS software. 

The spectra were fitted to a set of models (power-law, black- 
body, blackbody plus power-law, and two blackbodies, all 
modified by photoelectric absorption) using XSPEC version 
11.3.1. The single component models gave only marginally 
acceptable fits in most observations, while better results were 
obtained with the two component models. The best- fit param- 
eters of the latter models are reported in Table 2. 

In order to check for flux variability, we have also simul- 
taneously fit the five available PN spectra and the two MOS 
spectra for observation C with a double blackbody model with 
all parameters linked to the same value and a variable normal- 
ization factor From this analysis we can exclude significant 
(>3 (t) flux variations larger than ^20% among the different 
XMM-Newton observations. 

Since no significant differences in the spectral parameters 
are detected and the calibration of the PN instrument has 
proven to be very stable throughout the XMM-Newton mis- 
sion (see, e.g., Kirsch et al. 2005), a cumulative spectrum of 

^ Only the data of the two first observations have already been published 
(Lamb et al. 2002; Naze et al. 2004; Majid et al. 2004; McGarry et al. 2005; 
Nakagawa et al. 2007); we reanalyzed them using more recent calibration 
files, in order to consistently compare the results with those of the new obser- 
vations. 
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Fig. 1.— Cumulative PN spectrum of CXOUJ010043.1-721134 during 
observations B, D, E, and F. The best-fit double blackbody model is also 
shown together with the residuals and the separate contribution of the two 
blackbody components (dotted and dash-dotted lines). 

the PN data of observations B, D, E, and F^ has also been ex- 
tracted. The resulting net exposure time is 58 ks. Only the 
double blackbody model gives an acceptable fit to the cumu- 
lative spectrum with kTi ~ 0.3 keV and kTq, ^ 0.7 keV (see 
Fig. 1 and Table 2). The hydrogen column density is in good 
agreement with the average value of Nh = 5.9 x lO^'' cm"-^ 
expected towards this region of the SMC (Dickey & Lock- 
man 1990). The observed flux in the 2-10 keV energy range 
is 1.4 X 10"^' erg cm"^ s~\ corresponding to an unabsorbed 
luminosity of 6.1 x lO'"* erg s"^ for a distance of 60 kpc. The 
double blackbody model gave also the lowest values for 
the spectra of the single observations, but in those cases it 
was not the only one compatible with the data. In particu- 
lar, the power-law plus blackbody model, that is rejected with 
high confidence by the fit to the cumulative spectrum, gave 
acceptable fits to all the single spectra. 

In addition to the phenomenological models described 
above, we have also fitted the cumulative spectrum with the 
magnetar spectral model described in Rea et al. (2008). This 
model, originally proposed by Lyutikov & Gavriil (2006), is 
based on cyclotron resonant scattering of blackbody radia- 
tion in a twisted magnetosphere (Thompson et al. 2002). Al- 
though the resulting Xr (1.20 for 100 degrees of freedom) is 
slightly worse than for the double blackbody model, the fit is 
acceptable. The photoelectric absorption (Nh = 5 ± 1 x 10^" 
cm ) and blackbody temperature {kT = 0.32 ±0.08 keV) 
are consistent with the values derived from the double black- 
body fit (Nh = 6.3^^;^ x lO^^ cm'^ and kT = 0.30 ± 0.02 keV 
for the cooler blackbody). The best- fit values of the other 
spectral parameters are a resonant scattering optical depth of 
Tres = 1 .2 ± 0.2 and a particle velocity of I3t = 0.48 ±0.12. 
These parameters are in the same range as the ones observed 
in the other magnetar candidates (Rea et al. 2008). Although 
a direct information on the size of the emitting region cannot 
be derived in the current version of this model, an approxi- 
mate estimate gives a radius similar to the one of the cooler 
component in the double blackbody model. 

By inspecting the residuals from the best-fit models, we 
found no significant absorption or emission narrow-line fea- 
tures. We computed upper Umits on narrow lines' equivalent 
widths as a function of the assumed line energy and width 
This was done by adding Gaussian components to the 

* Observation A has been excluded because it was taken in a different 
operating mode, while no PN data were available for observation C. 
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TABLE 2 

Summary of the EPIC Spectral Results in the 0. 1-10 keV Energy Range. Errors Are Given at the 90% Confidence Level. 



Observation 


Model 


(lO^^cm-^) 


■n 
i 


PL norm.'^ 


kBlBBl 

(keV) 


Kbbi 
(km) 


I^bTbbi 
(keV) 


Rbbi!^ 
(km) 


Xr (d.O.t.) 


A 


PL+BB 




1 3+11 


0.8^:? 


"■--'^-0.03 








1.18 (83) 




BB1+BB2 




"■ -0 09 


12.9!^5^7 




3.6!i:i 


1.11 (83) 


B 


PL+BB 




1 7+0'' 
^ • ' -0.5 


2.o!f:^ 




o c+l.S 
°--'-0 6 


0.93 (115) 




BB1+BB2 








0.31 ±0.03 


11.0!f;« 


0.7±0.1 




0.72(115) 


C^ 


PL+BB 




~2 


<16 


0.34!J]:0^ 


9.6!|8 




1.16 (107) 




BB1+BB2 


-'■-'_4 4 






<0.39 


<11.6 


~1 


<3.7 


1.16(107) 


D 


PL+BB 


o ^+6.6 
J.O_35 


1 6+"'' 


2 0+5-0 


0.36 „ n3 


8 2+'0 






1.07 (83) 




BB1+BB2 


< 3.6 


^5+<)-02 

-!!-»3 


8 9+^^:^ 


0.9!0;1 


7+0-' 
"■'-0.4 


1.02 (83) 


E 


PL+BB 


4^2+7-2 

4,4^^:? 

7 6+^^^ 
"■"-3.4 




1 5+50 


37+fl-02 
-fl.Q3 




0.82 (84) 


F 


BB1+BB2 

PL+BB 
BB1+BB2 


1 6+20 


2.6±2.1 


33+0 0' 
-Q.Q6 

3^+<'-03 
-Q.Q5 

31+<'-0* 


io.i!?:5 

Q Q+2.6 

1 1 7+5-'! 




1 2+' -6 
^■^-0.7 

^■"-0.6 


0.75 (84) 
1.39 (65) 
1.30 (65) 


B+D+E+F* 


PL+BB 


9 1+7.9 


2 0+0-5 
^•"-0.4 


3 7+5-5 
^•'-1.9 


0.34 ±0.02 


O+0.9 
-0.,7, 


0.68+0:0^ 


1.75 (100) 




BB1+BB2 


6 3+^:^ 
"•■^-1.6 






0.30 ±0.02 


1^-1-1.4 


^•'-0.5 


1.14(100) 



^A free normalization factor has been introduced to account for inaccurate flux reconstruction in the MOS spectra, where, in most cases, the source is located 
on a CCD gap. 

''Assuming photoelectric absorption cross section from Balucinska-Church & McCammon (1992) and abundances from Anders & Grevesse (1989). 
''In units of 10"5 ph. cm"^ s"' keV"i, at 1 keV. 
"i Assuming a distance of 60 kpc. 
^Only MOS data. 
fQnly PN data. 




-loo r , 1 1 

0.5 1.0 1.5 2.0 2.5 

Line energy (keV) 

Fig. 2. — Upper limits (at 3 cr) to the equivalent width of narrow 
lines (either in emission or absorption) in the PN cumulative spectrum of 
CXOUJ010043.1-721134 during observations B, D, E, and F. 

double blackbody model and computing the allowed range in 
their normalization. The results for the high-statistics cumula- 
tive PN spectrum are summarized in Fig. 2, where the plotted 
curves represent the 3 <j upper limits for (j£ = eV. 

We performed a timing analysis to measure the source pulse 
period in each data set. After correcting the photon arrival 
times to the Solar system barycenter, we derived the best pe- 
riod values based on a periodogram analysis (Buccheri 
et al. 1983). The resulting values are indicated in Table 1. 
Considering also the periods measured by Chandra (McGarry 
et al. 2005), a linear fit to the ten values yields a period deriva- 
tive P = (1.9 ±0.1) X 10-11 g g-i (^2 of 1.32 for 8 degrees of 
freedom). 

Since observation D and E were performed only two days 
apart, we tried to better constrain the spin-down rate through a 
phase-coherent timing analysis of the two datasets. However, 
the periods uncertainties during each observation are too large 
to allow the prediction of the phase of the next observation to 



better than a pulse cycle. 

Many AXPs and SGRs are known to exhibit significant 
changes in their pulse profiles (e.g. Kaspi 2007; Gogiig et al. 
2002). To search for possible pulse shape variations in 
CXOUJ010043. 1-721 134 as a function of time, we com- 
pared the folded light curves using a Kolmogorov-Smimov 
test. Taking into account the unknown relative phase aUgn- 
ment, all the light curves are compatible with the same profile. 
We therefore summed them after appropriate phase shifts^. 
The resulting pulse profiles in the soft (0.2-1 keV, S) and hard 
(1-6 keV, H) energy ranges, together with their hardness ra- 
tio (computed as {H-S)/(H+S)), are shown in Fig. 3. This 
analysis does not show any significant profile changes with 
energy. The pulsed fraction^ in the 0.2-6 keV energy range is 
(32 ±3)%. 

3. DISCUSSION 

The new XMM-Newton observations reported here indicate 
that CXOUJ010043. 1-721 134 continued to spin-dow at a 
rate of ~1.9 x lO'i^ s s~^, consistent with previous results 
(McGarry et al. 2005). Although the errors on P are quite 
large, there is no evidence for major torque changes. Also the 
pulse profile, spectral shape and flux are consistent with no 
major changes, confirming that this AXP is one of the most 
stable members of its class. This characteristic allows us to 
sum up all the XMM-Newton observations taken with the same 
instrumental settings. 

The cumulative PN spectrum of CXOU JO 10043. 1-721 134 
cannot be adequately fit by a power-law plus blackbody 
model. This is the first magnetar for which such a conclu- 
sion can be derived based only on the X-ray data, thanks to 
the very low interstellar absorption. 

A good fit is instead obtained with a double blackbody 

5 We selected the shifts that maximized the Kolmogorov-Smimov statistics 
comparing subsequent observations. 

6 The pulsed fraction is defined as (Cmax-Cn,i„)/(Cmax+Cn,in), where C^ax 
and Cmiii are the background-subtracted count rates at the peak and at the 
minimum. 
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Fig. 3.— Background subtracted pulse profile of CXOU J010043.1- 
721 134 in different energy ranges (as indicated in the panels) and correspond- 
ing hardness ratio (PN data of observations B, D, E, F). 

model. The known distance of the SMC (60 kpc, Keller & 
Wood 2006) allows a precise measure of the size of the emit- 
ting region of the two blackbodies. The ~2 km radius of the 
region associated to the hotter component is compatible with 
a hot spot on the neutron star surface. The radius of the emit- 
ting region of the cooler blackbody, 12. l;!;^ ^ km, is consistent 
with a large fraction of the magnetar surface. However, the 
strong pulsation below 1 keV (see Fig. 3), where this com- 
ponent dominates (see Fig. 1), indicates that it cannot come 
from the whole neutron star surface. A similarly large black- 
body radius (^10 km) was also derived from the spectrum of 
the AXP XTE J1810-197 observed by ROSAT before the on- 
set of its outburst in 2003 (Halpem & Gotthelf 2005) and from 
phase-resolved spectroscopy of XMM-Newton observations of 
the same object (Israel et al. 2008) and of the AXP IE 1048. 1- 
5937 (Tiengo et al. 2005). However, in these cases, the less 
accurately known distance and the high interstellar absorption 
produce large uncertainties on the emitting region size. 

Assuming that the thermal photons are produced on the 
neutron star surface (and not, for instance, in the magne- 
tosphere) and considering that the blackbody is the most 
efficient thermal emission process at a given temperature, 
the radius of the region emitting the colder blackbody in 
CXOU JO 10043. 1-72 1134 is a firm lower limit to the radius 
of the compact object. This limit is not large enough to ex- 
clude any of the most popular equations of state for neutron 



stars, but it is the most constraining lower limit ever obtained 
for a magnetar. 

Magnetar spectra are expected to be more complex than 
a double blackbody. In fact most magnetar spectra cannot 
be fitted by such a simple model, that underestimates their 
emission in the 5-10 keV energy range. In all these cases a 
power-law tail in hard X-rays (>20 keV) has been detected 
(see, e.g., Kuiper et al. 2006; Gotz et al. 2006; Leyder et al. 
2008) and it is likely responsible also for the hard excess be- 
low 10 keV (Nakagawa et al. 2007). From a theoretical point 
of view, the emission expected from a magnetar has two com- 
ponents, with a thermal part directly from the surface, and a 
non-thermal one due to emission reprocessed in the magne- 
tosphere (Lyutikov & Gavriil 2006; Fernandez & Thompson 
2007). We found that also CXOU J010043. 1-721134 can be 
fit by a model of this kind (Rea et al. 2008). 

No compelling detections of spectral features in the per- 
sistent X-ray emission of magnetars have been reported so 
far. XMM-Newton and Chandra observations yielded strong 
upper limits (equivalent width <10 eV) for 4U0142H-61 
(Juett et al. 2002), IE 1048.1-5937 (Tiengo et al. 2005), and 
1E2259H-586 (Woods et al. 2004) in the 0.7-5 keV range. 
However, the high interstellar absorption towards these ob- 
jects causes a series of absorption edges in the observed spec- 
trum at low energies, introducing large systematic uncertain- 
ties in the search for features in the intrinsic AXP spectrum. A 
hint of a spectral feature at ^0.9 keV was noted by Durant & 
van Kerkwijk (2006) in the spectrum of the AXP 4U 142H-6 1 , 
after its deconvolution from interstellar absorption using the 
edges directly observed in the high resolution X-ray spectra, 
but also this measure is affected by the poorly constrained 
abundances of most interstellar elements. As shown in Figure 
2, for CXOU J010043. 1-721 134 we did not find evidence for 
lines, but, although this source is one of the dinmiest AXPs, 
we could put stringent hmits on narrow features in the soft 
X-ray band, which are virtually independent of the photo- 
electric absorption model. The dipolar magnetic field de- 
rived from the spin-down rate of CXOU JO 10043. 1-72 11 34 
is 4 X 10^^ G, corresponding to a proton cyclotron energy of 
~2.5 keV; however, a cyclotron line at lower energies is ex- 
pected if the cyclotron emission or absorption process occurs 
far from the neutron star surface, while a feature at higher 
energies is produced if strong multipolar magnetic field com- 
ponents are present. These effects, in addition to other pro- 
cesses that can suppress the spectral features (see, e.g.. Ho 
& Lai 2003), make the lack of proton cyclotron lines in the 
X-ray spectra of magnetars compatible with their magnetic 
fields of (6-250) x 10^^ G (corresponding to proton cyclotron 
energies of 0.4-15 keV) derived from their timing properties. 

We thank Nanda Rea for fitting the spectrum with the RCS 
model and useful discussion. This research has been partially 
supported by the Italian Space Agency. 
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